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Abstract 

We calculated the color-octet contribution to the J/tp hadroproduction at 
fixed target energies -\/T s ) — 40 GeV. We consider the J/i/j production with 
transverse momenta which can not be explained by primordial motion of par- 
tons, px > 1-5 GeV. It is shown that color octet contribution is dominant at 
these energies and reduces large dicrepancies between experimental data and 
color singlet model predictions. Having taken into account both contribution 
one needs a K-factor about 2-3 to explain the experimental data. 
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Investigation of the processes of heavy quorkonia production gives an excellent 
possibility to study a perturbative feature of QCD and allows one to extract long- 
distanse effects connected with hadronization phaze. Due to the large mass of c 
and b quarks production of the heavy quark-antiquark pair takes place at the short 
distances (~ tuq) and can be controled within the framework of perturbative QCD 
(PQCD). Then, the quark pair is bound into a quarkonium in a time scale of order 
of the the inverse binding energy, r ~ l/(m,QV 2 ), where v is the velocity of quarks 
within the bound state. As uiq — > oo, the heavy quark velocity inside the bound state 
v ~ l/ln(m q ) — > 0, and above mentioned two time scales become widely separated 
rriQ > rriQV 2 (for 66-system v 2 ~ 0.08 and for cc - v 2 ~ 0.23 — 0.3, charmonium is not 
truly nonrelativistic system). The separation of time scales makes possible to factorize 
production process of a heavy quarkonium in the short distance and long distance 
parts. Short distance part describes the heavy quark-antiquark pair production and 
can be controled perturbatively. 

The long distance part is related to matrix element responsible for the hadron 
formation from the quark-antiquark pair. 

Quarkonium production has traditionally been calculated in the color singlet 
model (CSM) |]. in this approach is proposed that quark-antiquark pair is produced 
in a color singlet state with the quantum numbers of the corresponding hadron. This 
heavy pairs eventually creates the hadronic state with a probability determined by 
the appropriate quarkonium wave function at the origin. It is assumed that for heavy 
quarks soft gluon emmision as well as other nonperturbative effects such as higher 
twist contributions are negligible. While this model gives a reasonable description 
of J /if) production cross section shape over p-p or if it completely fails in the ex- 
planation of the integrated cross cross section (K factor 7-10 is needed to explain 
experimental data). The anomalously large cross section of J/if> production at 
large transverse momenta at the Tevatron reveals another bad feature of the CSM. 
Within the framework of the CSM it is impossible to explain the anomalously large 
ip H and direct J /if) production [Q at CDF experiment. 

Indeed, the requirements of the CSM is very strong and lead to the suppression 
of both fusion and fragmentation contributions. Production of a quark-antiquark 
pair on the time scale r ~ I/tuq with the proper quantum numbers of subsequent 
hadronic state must be accompanied with the emmition of hard gluons. Therefore, 
production cross section in the CSM is suppressed by powers of a s /Ti. 

The CSM is a nonrelativistic model where the relative velocity between the heavy 
constituents in the bound state is neglected. But discrepancies between experimental 
data and CSM predictions hint that 0(v) corrections as well as other mechanisms of 
quarkonia production, which do not appear at leading order in v, should be considered. 
Treating higher order in v it is possible to avoid the supression caused by emission 
of hard gluons, if one allows a quark-antiquark pair to be produced on the time scale 
1/rriQ in any color state (singlet or octet) with any quantum numbers. Such a pair 
evolves into a hadron state by emitting of soft gluons with momenta of an order of 
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itiqv 2 . Expansion of quarkonium cross sections and decay widths in the powers of 
relative velocity v of heavy quarks in a bound state has recently been realized in 
terms of Nonrelativistic QCD (NRQCD) |J. This formalism implies not only color- 
singlet processes but so called color-octet mechanism, when a quark- ant iquark pair 
is produced on the small time scales in the color octet state and evolves into hadron 
by emmition of soft gluons. According to the factorization approach based on the 
NRQCD, the production cross section for quarkonium state H in the process 

A + B -> H + X (1) 

can be written as 

Gi i = Yl J dx 1 dx 2 fi/A(x 1 )f j /B(x 2 )cr(ij -> H) (2) 
id 

a(ij -> H) = J2C t] [n] < 0\O H [n]\0 > 



where the fi/A is the distribution function of the parton i in the hadron A. Sub- 
prosess cross section is separeted into two: short distance, (7 ,J [n], and long distance, 
< 0|O^[n]|0 >, parts. The C y "[n] is the production cross section of a heavy quark- 
antiquark pair in the i and j parton fusion. It should be calculated in the framework of 
PQCD. The [n] state can be either a color sing let or octet state. The < 0\O H [n}\0 > 
describes evolution of a quark-antiquark pair into a hadronic state. These matrix 
elements cannot be computed perturbatively. But relative importance of long dis- 
tance matrix elements in powers of velocity v can be estimated by using the NRQCD 
velocity scaling rules |J. 

The shapes of the px distribution of short distance matrix elements within the 
color octet model indicates that new mechanism can explain the Tevatron data of 
direct J/ip and ip production at large px- But unlike color-singlet matrix elements 
connected to the subsequent hadronic nonrelativistic wave functions at the origin, the 
color octet long distance matrix elements are unknown and should be extracted from 
experimental data. Explanation of the CDF data of S state charmonia production 
at the large px will be succesfull after testing the values of color octet long distance 
matrix elements in other experimental data 

The color octet contribution to the J/ip photoproduction was analized in the pa- 
pers |7|, ||. Recently, the J/tj) hadroproduction at fixed target energies has been 
studied by including color-octet mechanism || [H|. Large discrepancies between ex- 
perimental data and CSM predictions for the J/ip production total cross section were 
explained. The color octet contribution is dominant in the J/ip hadroproduction at 

energies y[s) ~ 30 — 60 GeV. The analises carried out in these papers J7| ||, 0, [10f| 
demonstrate that to fit the photoproduction and hadroproduction data at the low en- 
ergies require rather small values for color octet matrix elements than those extracted 



from a prompt J / ip production at CDF [12 
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In the present letter we consider the color octet contribution to the J / ip hadropro- 
duction with nonzero px at fixed target energies. We called nonzero px such trans- 
verse momenta which cannot be explained by primordial motion of partons inside 
the colliding hadrons. We consider px > 1.5 GeV to calculate the short distance 
matrix elements perturbatively and avoid the infrared and collinear divergencies in 
the production of 1 Sq and 3 Po,2 quark-antiquark states. 

In the papers || [H| for calculation of the total cross section of the J/xjj produc- 



tion only the subprocesses 2 — ► 1 were taken into account. These subprocesses are 
the lowest order in perturbative series over a s and give the main contribution to the 
integrated cross section. The transverse momenta of produced particles due to inter- 
nal motion of partons are of an order A.qcd, and pt behaviour of the cross section 
cannot be controlled perturbativle. The subprosses 2 — > 2 are the lowest order which 
contribute to the J /if) production at px » AqCD. 

99 -> (cc)g, 

gq -> (cc)q, (3) 
qq -> (cc)g, 

where (cc) denotes any state of heavy quark-antiquark pair. These subprocesses give 
small contribution to the total cross section because additional a s and the steeply 
falling down px behaviour. The cross section of the J/ip production can be written 

as 

°Jli> = v(J/^)dir + E Br (Xc.j -> J/^X)a XcJ + Br(ip' -> J/^X)a^, (4) 

J=0,l,2 

where in the production of each quarkonium state is contributed by both color singlet 
and octet states, 

a( J^) dir = a ^ + f m + ^ v(QQ( 2s+1 P!j 8) ) < 0\O J ^C^Lj\0 > (5) 

where the sum stands over the states 3 -Po 12; 1 Sq and 3 5*^. The expressions for dif- 
ferential cross sections for color singlet states we take from 11 1. For color octet 



states the short distance matrix elements have recently been calculated by Cho and 



Leibovich 12 



As concerned the long distance matrix elements their number should be reduced 
by using NRQCD spin symmetry relations: 

< 0|Of ( 3 Pj)\0 >= (2J+ 1) < 0|Of ( 3 P )|0 >, (6) 
< 0\Of J ( 3 Si)|0 >= (2J+ 1) < 0|Of°( 3 Si)|0 > . (7) 

After that from the parameters which give main contribution in the cross section 
only four independent matrix elements remain; < 0|Og^( 3 S'i)|0 >, < 0|Og cl ( 3 S'i)|0 >, 
< 0\O^{ 3 P )\0 > and < 0|Og /,/ '( 1 So)|0 >• Unfortunatly, values for these matrix 
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elements obtained from fitting various experimental data are different. There are 
two different values for the matrix elements < OlOg^C 3 ^) |0 >, < 0|Of 1 ( 3 S 1 )\0 >, 
extracted from CDF data by two group of authors. Using only the dominant frag- 
mentation contributions to the J/ip production Caciari et al.Q obtained the following 
values: 

< OlOg^fSOlO >= 15 ■ 10- 3 GeV 3 , 

< 0|Of 1 { 3 S 1 )\Q >= 2.4 • 10~ 2 GeV 3 . 

Using the full perturbative expressions for the short distance matrix elements Cho 
and Leibovich obtained smaller values by factor two for the above parameters, 6.6 • 
10- 3 GeV 3 and 9.8 • 1(T 3 GW 3 , respectively. 

As for the other two parameters it is possible to extract only their combinations 
from experimental data. From charmonium production at large transverse momenta 



at CDF Cho and Leibovich extracted [112 
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< 0\O J 8 /lp ( l S )\0 > +^2 < 0|O 8 W ( 3 P )|0 >= 6.6 • 10- 2 GeV 3 . 



Different combination of these parameters were extracted from J / ip photoproduction 
data and hadroprodaction data at fixed target energies: 

< 0|0 8 J/,/ '( 1 So)|0 > +^ < 0|0 8 J/,/, ( 3 Po)|0 >= 2 • l(T 2 GeU 3 [10], 



< 0|O 8 J/ YSc)|0 > +^ < 0|O 8 W ( 3 P )|0 >= 3 • 10- 2 GeU 3 B |. (9) 



m 2 c 



As one can see from equations (8) and (9) the photoproduction and hadroproduction 
data at low energies are consistent with each other within the model error. But there 
is the large discrepancy between the low, (9), and high energy, (8), values of the 
parameters. Solving the system of equations (H) and @, one will obtain a negative 
value for the < 0^( 3 Pq) > matrix element. The reason for such discrepancy may 
be an overestimation of the value obtained in the paper [T2j (equation ||). In the 
prompt J/ip production at the Tevatron the 1 S'q 8 ' ) and 3 Pj 8 ^ states give a dominant 
contribution at transverse momenta near px = 5 GeV. The value 6.6 ■ !0~ 2 GeV 3 was 
extracted just from this region of pt- But at these values of transverse momenta J ftp 
is produced mainly at small partonic x (x ~ 1CT 2 -j- ICT 3 ) in the gluon-gluon fusion. 
Cho and Leibovich in paper JL^] use the partonic parametrization MRSDO ||14|| . But 



in the above mentioned region of x the MRSDO parametrization gives rather smaller 
value for the gluon distribution function than more realistic parametrization (GRV 
LO or GRV HO Using GRV LO parametrization gives 1.55 times large values 



for the production cross sections of 1 S and 3 Pj states. So, more realistic value for 
the combination is 

< 0|O 8 W ( 1 5 )|0 > < O|0^( 3 Po)|O >= 4-r 4.4 • 10- 2 GeV 3 . (10) 

mi 
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Consequently, value of 1.4 4 1.5 times large is needed for the < Og^( 3 Si) > matrix 



element to explain the prompt J/ip production cross section at px — 10 GeV ||12| , 

< Oi'^^Sx) >= 9 4 10 ■ 10- 3 GeV 3 . (11) 

After these changes the dicrepancy between the combinations (9) and (10) would not 
be large for the radical choice < Oq^( 3 Pq) >= 0. 

We presented here, fig.l, results for the differential cross section for the set of 
parameters taken from [fl^| . We asumed that < O^^Sq) >=< 0$^( 3 Po) >. 

In the fig. 2 presented the differential cross section for another set of parameters 
obtained after taking into account above mentioned corrections: 

< Og^fSO > = 10.5 • 10- 3 GeV 3 , 
<0\O% cl ( 3 S 1 )\Q> = 9.8 • 10- 3 GeV 3 , 
<O J 8 /lp CS )> = 3.7-10- 2 GeV 3 , 



<o 8 "T^o)> = o. 



(12) 



For < Oq^^Sq) > we choose the medium value of the combinations (|S[) and (|TUD . For 
parameters which describe the transition color octet state into a -?//-meson we use the 
values from Theoretical predictions are compared with the experimental data 
of experiment E689 at FNAL |13] . We use the functional form of for px distribution 
from |l3j to reproduce the curve of experimental data (fig. 1,2). As one can see from 



figures color octet contribution is dominant. But the total theoretical prediction 
are small and K-factor about 2 4 2.5 is needed to explain the experimental data. 
It is wort to mention that in our calculations there is uncertainty that lead to the 
decreasing of the theoretical prediction in the whole region of considered pt- The 
relative velocity of c and c quarks in the charmonium is about v 2 ~ 0.23 4 0.3. 
This means that soft gluons emmited from color octet quark- ant iquark state while 
evolution into J/ip has an momentum about 0.7 4- 1 GeV (2m c v 2 ). So, it is necessary 
to produce color octet state of quark-antoquark pair with mass large then 2m c . In 
the case of the photoproduction or the total cross section in hadroproduction at fixed 
target energies the kimematical effect from the difference between the mass of J/if> 



and color octet quark- ant iquark pair is very large [JTUj] since the gluon distribution 
rises steeply at small x. This reduces cross sections twice and 'true' matrix elements 
would therefore be large than those extracted by using the small mass of quark- 
antiquark pair fl0| . In our case the influence of such effect is not so large and leads 
to the decreasing the cross section by about 25% at p t ~ 1.5 GeV and by about 10% 
at pt ~ 3 GeV. Another uncertainty comes from decay process of quark-antiquark 
pair into J /if}. The J /if} production cross section shape over px differs from that 
of color-octet quark-antiquark pair. For qualitative estimations of this corrections 
we assumed that produced quark-antiquark state is unpolarized and its decay into 
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J/ip and gluon is isotropic. In this approach we calculated the deviation of the J/ip 
Pt distribution shape from the distribution shape of heavy quark- ant iquark pairs 
with mass 4 GeV. The difference is negligible near p T ~ 1.5 GeV and at 3 GeV is 
about 20%. Unfortunately, we can only qualitatively estimate these two corrections 
because the true mass of a quark- ant iquark state before evolution into hadronic state 
is unknown. But it is obvious that these two corrections lead to the reduction not to 
the rising of the cross section. 

In conclusion, we have calculated the color octet contribution to the J/ ip hadropro- 
duction with nonzero pr at fixed target energies. Color octet processes are dominant 
in the J/ip production and large about an order of magnitude then color singlet 
contribution. 

After taking into account color octet contribution the K-factor about 2 ~ 3 is 
need to explain the experimental data. Uncontrolled corrections coming from the 
decay process of quar- ant iquark states into J/ip and 'soft' gluons leads only to the 
decreasing of the cross section and discrepancies become larger. 

We are indebted to N. Kochelev, W.-D. Nowak and O. Teryaev for useful discus- 
sions and helpful comments. 
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FIGURE CAPTIONS 



Fig.l. Transverse momentum diferential cross sections for the set of color octet matrix 
elements from ||12|| ; solid curve represents experimental data from [|13|], dashed curve 
- total theoretical predictions and dotted curve - CSM predictions. 

Fig. 2. Transverse momentum diferential cross sections for the color octet matrix 
elements from eq.(12); The curves in this figure are labeled in the same way as in 
fig.l. 
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